CERN-TH/96-285 



hep-ph/9610350| 



HIGLU AND HDECAY: Programs for 
HiGGS Boson Production at the LHC and 
HiGGS Boson Decay Widths* 

Michael Spira 

CERN, TH Division, CH-1211 Geneva 23, Switzerland 



Abstract 

The total cross section for Higgs boson production via the dominant gluon fu- 
sion mechanism including NLO [two-loop] QCD corrections can be calculated 
numerically with the program HIGLU. The QCD corrections are included for 
arbitrary Higgs and quark masses and increase the cross section at the LHC 
by up to a factor of 2. The source code HIGLU provides the evaluation of the 
production of the Standard Model [SM] Higgs boson as well as the neutral 
Higgs bosons of the minimal supersymmetric extension [A4SSAi]. 

The program HDECAY determines the decay widths and branching ratios 
of the Higgs bosons within the SM and the AiSSM, including the dominant 
higher-order corrections. The latter are dominated by QCD corrections and 
two-loop corrections to the couplings and Higgs masses of the A4SSA4. The 
program includes all decay modes with branching ratios larger than 10~*. 
Moreover, below-threshold decays with off-shell top quarks, gauge and Higgs 
bosons are implemented. In addition the program is able to calculate the 
branching ratios of the MSSM Higgs bosons into supersymmetric particles, 
which can be dominant. 
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The total cross section for Higgs boson production via the dominant gluon fusion mech- 
anism including NLO [two-loop] QCD corrections can be calculated numerically with 
the program HIGLU. The QCD corrections are included for arbitrary Higgs and quark 
masses and increase the cross section at the LHC by up to a factor of 2. The source 
code HIGLU provides the evaluation of the production of the Standard Model [SAI] 
Higgs boson as well as the neutral Higgs bosons of the minimal supersymmetric exten- 
sion [MSSM]. 

The program HDECAY determines the decay widths and branching ratios of the Higgs 
bosons within the SM and the MSSM, including the dominant higher-order correc- 
tions. The latter are dominated by QCD corrections and two-loop corrections to the 
couplings and Higgs masses of the MSSM. The program includes all decay modes with 
branching ratios larger than 10~*. Moreover, below-threshold decays with off-shell top 
quarks, gauge and Higgs bosons are implemented. In addition the program is able to cal- 
culate the branching ratios of the MSSM Higgs bosons into supersymmetric particles, 
which can be dominant. 

1. Introduction 

At the future CERN hadron collider, the LHC, the Higgs boson will be produced 
primarily via the gluon fusion mechanism for the entire relevant Higgs mass range 
within the Standard Model [iSA^] u as well as its minimal supersymmetric extension 
[AliSiSA^] □. The two-loop QCD corrections to the production cross sections of 
scalar [CP-even] as well as pseudoscalar [CP-odd] Higgs bosons have been calculated 
in uu for heavy squarks. The QCD corrections are significant for the theoretical 
prediction of the cross sections leading to an increase by up to a factor of 2 compared 
with the lowest-order results. Recently, the leading part of the pure QCD corrections 
to the squark loop contributions to MSSAA Higgs boson production have also been 
evaluated by means of low-energy theorems in the limit of heavy gluinos □. The 
squark loops are sizeable for squark masses below ^ 400 GeV. The corrections are 
of the same magnitude as the corrections to the quark loops so that the K factors 
are only slightly changed by the inclusion of the squark loops. 

The program HIGLU^provides the calculation of the total Higgs production cross 
sections including next-to-leading order QCD corrections to the heavy quark loops I. 
Various input parameters can be chosen from an input file, including a flag specifying 
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the model. Possible options are the Standard Model, its minimal supersymmetric 
extension, and a general Higgs model by initializing the Higgs Yukawa couplings 
appropriately. The program includes the contribution of the top and bottom quarks 
in the loops that generate the Higgs couplings to gluons. Within the Standard 
Model, as well as in most of the parameter space of the A4SSA4, these contributions 
provide an excellent approximation for all practical cases. Moreover, the program 
allows the calculation of the decay widths of Higgs bosons into gluons including 
next-to-leading order QCD corrections. The giuonic decay mode plays a significant 
role in the intermediate mass range at future e~^e~ colliders au. 

The search for Higgs bosons at the LHC mainly proceeds via looking for ZZ, 
W^W^ and 77 final states as decay products of the SM Higgs particle and in 
addition for Higgs boson pairs in the case of the AiSSAd Higgs bosons. Thus it is 
of vital importance to have reliable predictions for the branching ratios of the Higgs 
bosons in both models. These are determined by the program HDECAY including 
all relevant higher-order corrections, which are dominated by QCD corrections and 
the two-loop corrections to the JASSM Higgs masses and couplings arising from 
the large top mass. HDECAY □ calculates all decay modes with branching ratio 
larger than 10~* and includes below-threshold decays involving off-shell top quarks, 
gauge and Higgs bosons in the final states U. Moreover, all decay modes of the 
A4SSA4 Higgs bosons into SUSy particles are contained in the program, which 
can be dominant for Higgs masses below the tt-threshold. 

The source codes of the programs are written in FORTRAN. They have been 
tested on computers running under different operating systems. The numerical 
ase space integration of HIGLU is performed by using the VEGAS package 
for integrals of dimension up to 3. Parton distributions can be attached to 
HIGLU in any desirable way by adjusting the corresponding subroutine. As stan- 
dard parametrization, the program contains the GRV sets til. 

2. Processes 



2.1. pp-^^^ + X 

The hadron cross section of Higgs boson'' production via gluon fusion gg 
(4>° = H, A) including [two-loop] QCD corrections, can be cast into the form 

^(pp ^ cE,o + X) = ^f^ + Aa*:* + Act*; + Aa*; + Aaf^ . (1) 

The lowest-order cross sections ct^q can be found in RefsE^i. The term Acr*"^ 
parametrizes the infrared regularized virtual two-loop corrections and the terms 
Act* = g,q,q) the individual coUinear regularized real one-loop corrections 

''The scalar [CP-even] Higgs particles will generically be denoted by the pseudoscalar [CP-odd] 
by A and all the neutral Higgs bosons by The A^iSiSA^ Higgs particles, including the charged 
ones, will be denoted by 
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corresponding to the subprocesses 

g5->$V gq^<^\. qq^^"g- (2) 

Their expressions can be found in Refs. 

The calculation of the [two-loop] virtual corrections requires the evaluation of 
massive two-loop three-point functions depending on one free parameter, the ratio 
of the Higgs and heavy quark masses. The five-dimensional Feynman integrals 
have been analytically reduced to one-dimensional integrals; they are numerically 
integrated by an adaptive Romberg integration 0. The regularization of UV, IR and 
coUinear divergences is performed in n = 4— 2e dimensions. The UV singularities are 
absorbed in the strong coupling and the on-shell quark mass mg. The infrared 
singularities of the virtual corrections cancel against the corresponding singularities 
of the real corrections of the subprocesses given in Eq. (||). The left-over collinear 
poles are mapped into the renormalized NLO parton densities. The renormalization 
and factorization schemes have been chosen as the MS scheme. 

The program HIGLU calculates the five terms in Eq. (|l|) contributing to the 
total cross section separately, as well as their sum, for all kinds of neutral Higgs 
bosons <i>°. 

2.2. 

The decay widths of Higgs bosons into gluons up to next-to-leading order 
are given by 



r($"-^5.g(5),5gg) = TLoi^^^gg) 



TT 



E^o = E'^,,,+E%^+NfE%^. (3) 

The coefficient denotes the infrared-regularized virtual two-loop corrections, 

E^gg and E^^g the infrared and collinear regularized real one-loop corrections. The 
analytical formulae of these contributiopA^pan be obtained from the gluon fusion 



process by crossing and are given in RefsJS'ollj. The parameter Np fixes the number 
of light external fiavors produced in the decay gqq, which has to be identified 

with the number of flavors contributing to the QCD f3 function in order to map 
large logarithms into the running strong coupling as{fi^). This issue is important, 
if bottom and charm quarks will not be included as external flavors in the gluonic 
decay mode, but added to the bb and cc decay channels: in order to avoid large 
logarithms in the perturbative expansion, the strong coupling has to be evolved with 
three light flavors instead of five at scales above the heavy-quark thresholds. The 
program HIGLU evaluates this decay mode, including the full mass dependence. 
The program HDECAY on the other hand only contains the leading part in the 
heavy-quark mass limit, which is a reliable approximation in all relevant cases. The 
gluonic branching ratios reach values of about 10-20% in the intermediate mass 
ranges and are thus significant for phenomenological analyses. 
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2.3. $" ^ QQ, H+ UD 

The expressions for decay widths of Higgs bosons into heavy-quark pairs can 
be found in Ref. IIJ. In order to absorb large logarithms in the limit A/$ 3> tiq, 
one has to introduce the running MS mass for the heavy quarks in the final states, 
which has to be evaluated at the scale of the corresponding Higgs boson mass M^. 
This definition leads to a well-behaved perturbative expansion in the case of charm 
and bottom quarks above the thresholds. The running MS mass is calculated 
from the value of the pole mass by means of the relation given in Rcf.E3. A flag 
is provided for the choice of the NLO or NNLO formula of this relation, which is 
important in particular for the charm-quark mass, because of a bad convergence of 
the perturbative expansion at scales of the order of the charm mass. The MS mass 
characterizes the preferred input parameter, because it can be obtained directly 
from fits to QCD sum rules Ej. 

In the limit M$ ^ mg the QCD corrections are known up to NNLO. At the 
threshold, mass effects are important. Thus we have included the full massive 
NLO QCD corrections in the program HDECAY, which are smoothly interpolated 
with the NNLO expressions in the large-Higgs mass range. The QCD corrections 
reduce the decay widths significantly. Electroweak corrections are small EZI and 
hence neglected in the program. The hh decay mode dominates in the intermediate 
neutral Higgs mass range, if decays into SUSy particles are kinematically forbidden, 
except for the heavy scalar MSSAA Higgs particle H for small tg/3. The decays 
into cc and t+t~ pairs can reach branching ratios of about 10%. Above threshold 
the tt decay modes are dominating in the MSSM. Charged MSSM Higgs bosons 
predominantly decay into tv and th pairs, if SUSy particle decays are impossible. 
For top quarks in the final states, HDECAY also takes into account below-threshold 
decays into off-shell top quarks cl, which are smoothly interpolated with the on-shell 
regime. 



2.4. $0 77 

The photonic decay widths of the neutral Higgs particles are built-up by charged- 
particle loop contributions from fermions and W bosons in the and in addition 
from charged Higgs particles, charginos and sfcrmions in the M.SSA/L The ex- 
pressions can be found in Ref.El. The QCD corrections are moderate u and thus 
neglected in the program HDECAY. The photonic branching ratios reach values of 
about 10^'^. These decay modes are important for the Higgs search at the LHC 0. 



2.5. n WW,ZZ 

The expresaons for the scalar decay widths into W and Z boson pairs can be 
found in Ref.E^I. They are the dominant decay modes of the SJ^ Higgs boson 
for Higgs masses Mh > 140 GeV. In the M.SSM. they are only important for 
small tg/3. The program HDECAY includes below-threshold decays into off-shell W 
and Z bosons which are interpolated with the expressions of the on-shell range. 
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Electroweak corrections are small L3 in the relevant Higgs mass ranges and neglected 
in HDECAY. 



2.6. Decays into Higgs particles 

The decaY_widths of the AiSSA4 Higgs bosons into lighter Higgs particles are 
given in Ref.Ej. The decay mode H ^ hh is dominant below the ti-threshold for 
small tg/3. In parts of the mass ranges, the decay modes A Zh and W^h 
are dominating. HDECAY also includes below-threshold decays into off-shell Z and 
Higgs bosons □ and interpolates these contributions with the on-shell regime. The 
dominant corrections are contained in the two-loop-corrected AiSSAd couplings. 

2.7. Decays into SUSy particles 

The partial decay widths into neutralinos, charginos and sfermions are given 
in Ref.E3. Their branching ratios can amount up to ^out 100% below the top 
threshold for small tg/3 and remain significant above E3. HDECAY evaluates all 
these partial decay widths within the JviSSAA^ in which they are completely fixed 
by two M.SSM parameters, /i and M2. 

3. Input and Ouput 

All masses, couplings and M.SSM parameters can be defined in separate input 
files for both programs HIGLU and HDECAY. Moreover, there are a few fiags to 
allow a choice of the kind of process, and of the conventions for the parameters 
involved in the different processes. Thus the two programs provide a very flexible 
and convenient usage, fitting to all options of phenomenological relevance. 

The program HIGLU writes all chosen input parameters as well as the results 
obtained by VEGAS integrations including the individual parts of the QCD correc- 
tions to the file unit 99. 

HDECAY produces many output files containing the results for the Higgs branch- 
ing ratios and total decay widths in table format. These can easily be read by other 
FORTRAN programs, which need the values for the Higgs branching ratios and 
decay widths. 
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